arXiv:1501.01191vl [nucl-th] 6 Jan 2015 


1 S'o nucleon-nucleon scattering in the modified Weinberg approach 

E. Epelbaum , 1 A. M. Gasparyan , 1,2 J. Gegclia , 1,3 and H. Krebs 1 

1 Institut fur Theoretische Physik II, 

Fakultdt fur Physik und Astronomie, 

Ruhr-Universitat Bochum 44 780 Bochum, Germany 
2 SSC RF ITEP, Bolshaya Cheremushkinskaya 25, 117218 Moscow, Russia 
3 Tbilisi State University, 0186 Tbilisi, Georgia 
(Dated: 12 July, 2014) 

Abstract 

Nucleon-nucleon scattering in the 1 S'o partial wave is considered in chiral effective field theory 
within the renormalizable formulation of Ref. (l|] beyond the leading-order approximation. By 
applying subtractive renormalization, the subleading contact interaction in this channel is taken 
into account non-perturbatively. For a proper choice of renormalization conditions, the predicted 
energy dependence of the phase shift and the coefficients in the effective range expansion are found 
to be in a good agreement with the results of the Nijmegen partial wave analysis. 
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I. INTRODUCTION 


The seminal work of Weinberg [2| has triggered a renewed interest to the nuclear force 
problem in the framework of effective field theory (EFT). In this approach, nuclear forces 
are defined as kernels of the corresponding dynamical equations and can be derived order- 
by-order making use of the systematic chiral expansion. 

Starting from the pioneering work of Ref. [3], this approach has developed rapidly and 
is nowadays commonly employed in studies of low-energy few- and many-nucleon dynamics 
and nuclear structure calculations, see jU-jfl for recent review articles. While offering many 
attractive features, Weinberg’s approach was criticized for being non-renormalizable. The 
main difficulty is related to the fact that iterations of the truncated NN potential within 
the Lippmann-Schwinger (LS) equation generate ultraviolet (UV) divergencies which cannot 
be absorbed by counter terms (contact interactions) included in the truncated potential. In 
particular, infinitely many counter terms are needed to absorb UV divergences emerging from 
iterations of the leading-order (LO) one-pion exchange (OPE) potential [t|. This feature is 
sometimes referred to as inconsistency of Weinberg’s approach. 

The UV cutoff A can be removed from the LS equation by enforcing the limit A —y oo non- 
perturbatively, see e.g. [S, §). It is possible to obtain a finite, manifestly non-perturbative 
solution of the LS equation with a singular l/r 3 -potential by including one/no contact 
operator in each attractive/repulsive channel [9j. However, such a procedure is incompatible 
with the principles of EFT which require that all UV divergences emerging from iterations 
of the LS equation are absorbed by counter terms [l3|. It is not surprising that such an 
approach fails to reproduce experimental data even at N 3 LO fill ]. 

Treating the exchange of pions perturbatively as suggested by Kaplan, Savage and Wise 
(KSW) jl2|] allows one to avoid the above-mentioned inconsistency. However, the perturba¬ 
tive series fail to converge within this framework fl3-16l . 

Presently, there exist different views and formulations of organizing the chiral expansion 
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32fl. A novel approach to the NN scattering problem in 


1, 3^-36] and is referred to as the modified Weinberg ap 


in the few-nucleon sector [8| 

EFT has been formulated in Refs, 
proach. Within this framework, the leading order (LO) NN scattering amplitude is obtained 
by solving the Kadyshevsky equation [37j] for the LO potential consisting of the contact inter¬ 
action part and the OPE potential. This equation provides an example of three-dimensional 
integral equations which satisfy relativistic elastic unitarity. An important feature of the 
Kadyshevsky equation is that it is renormalizable for the LO potential, i.e. all ultraviolet 
divergences generated by iterations can be explicitly absorbed into redefinition of the NN 
derivative-less contact interaction. The scattering amplitude can still be renormalized if 
higher-order corrections to the potential are taken into account perturbatively. If higher 
order corrections to the potential indeed provide small contributions to the amplitude, their 
perturbative and non-perturbative inclusions are expected to lead to small differences in 
the results which are beyond the accuracy one is working at. However, this observation 
is only meaningful if a proper renormalization is carried out in both cases. In general, we 
are not able to subtract all divergences from amplitudes if higher-order contributions in the 
potential are treated non-perturbatively. In the 1 S'o partial wave, one observes a very large 
discrepancy between the LO EFT results and the experimental data already at rather low 


energies [l[. This large discrepancy signals that at least a part of the higher-order contribu¬ 
tions in the effective potential is likely to require a non-perturbative treatment within our 
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approach. 1 

In this paper we study in detail the role of the next-to-leading order (NLO) short- 
range contribution to the potential which can be included both perturbatively and non- 
perturbatively. Specifically, we will express the solution to the integral equation in a closed 
form following the lines of Ref. (doj and apply the BPHZ-type subtractive renormalization 
[4lj ] . After subtracting all ultraviolet divergences, we will calculate the remaining finite 
expressions numerically, fit the available two low-energy constants (LECs) to the data and 
compare the obtained results with the phase shifts for various choices of the renormalization 
scale parameter. Here and in what follows, the resulting NN amplitude will be referred to 
as NLO as opposed to the LO result of Ref. [lj]. A more complete calculation including the 
corresponding two-pion exchange potential to Erst order in perturbation theory is postponed 
for a future study. 

Our paper is organized as follows: In section |TT] we provide the formal expression for 
the scattering amplitude by making use of the standard two-potential formalism. Subtrac¬ 
tive renormalization of the amplitude is discussed in detail in section UTT1 Next, section HVl 
addresses the issue of the appropriate choice of the renormalization conditions (i.e. subtrac¬ 
tion scale) and also presents the results of our calculation. Our Endings are summarized in 
section El 


II. FORMAL EXPRESSION FOR THE SCATTERING AMPLITUDE 


In the framework of the modified Weinberg approach, the NLO 1 Sq partial wave NN 
scattering amplitude 2 can be obtained by extracting the S-wave component from the solution 
to the integral equation (for the fully ofi-shcll amplitude T) 


T(p 0 ,p',P ) = V(p',p) + J d 3 kV (p',k)c{p 0 ,k)T(p 0 Xp) , (2.1) 

777 2 1 

G(po,k) = ^ -r. ( 2 - 2 ) 

2 ( 2 tt ) (p 2 _|_ m 2 \ _ y/k 2 -m 2 + ie) 

where p (p r ) is the incoming (outgoing) three-momentum of the nucleon in the center-of-mass 
frame, po = \J q 2 + m 2 with m denoting the nucleon mass and q being the corresponding 
three-momentum of an incoming (on-mass-shcll) nucleon. Further, the potential is given by 


V (p',p) 


C 


[C + C 2 (p' 2 +p 2 )] 
Vc + K, 


g\M 2 

4R2 ip, 


Cs — 3 Ct + 



+ D M 2 . 


1 

pf + Ml 


(2.3) 


Here qa , F n and M n are the nucleon axial-vector coupling, pion decay constant and the pion 
mass, respectively. The parameters Cs, Ct, C 2 and D refer to the LECs of the efiective 


1 Notice, however, that the LO calculations reported in Refs. 38, |3!J] within the standard nonrelativistic 


approach using a coordinate-space regularization for the OPE potential yields a superior description of 
the phase shift. 

2 Note the different overall sign in comparison with the Feynman amplitude considered in Ref. |lj]. 
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Lagrangian. Below, we work with the S'-wave component of Eq. (Eli and denote the 1 So 
partial wave projected OPE potential by 14 (p',p) with 


Vn(p',p) 


g\Ml (p - p'Y + M 2 
16 F 2 pp' 11 (p + p') 2 + M 2 ‘ 


(2.4) 


For the analysis of divergent integrals, it is useful to have the asymptotics of V 7T (p',p) at 
large values of momenta 


vAp',p)\ 


p—t oo, p < oo 


9a m I 

4 F 2 p 2 


K(p + 1,p) I 


>oo, |ZI<oo 


g\Ml 

8 


In p. 


(2.5) 


The contact-interaction part of the potential Vc is separable. Therefore, it is possible 
to write the solution to Eq. (12. ip in a form, which allows one to carry out the subtractive 
renormalization explicitly. This can be achieved by making use of the well-known two- 
potential formalism. For this purpose, we write Eq. (12.ip symbolically as 


T = V + V GT, 

(2.6) 

and express its solution as 


T = T n + (l + T n G)T c (l + GT n ), 

(2.7) 

where T n and Tq satisfy the equations 


Tit = 14 + 14 G T n , 

(2.8) 

T c = V c + V c G(l + T n G)T c . 

(2.9) 

For a separable contact-interaction potential, 


Vc(p',p ) = £(p’) T C£(p), 

(2.10) 


where C and £(p) are 2x2 and 2x1 matrices, respectively, whose explicit form will be 
specified below, the solution to Eq. (12.9p is also given in a separable form 

Tc(p 0 ,p',p)=f(p')xap), (2.11) 

where X is a 2 x 2 matrix, 

X = [C -1 — E] -1 , (2.12) 

and the 2x2 “selfenergy” matrix E reads 

£(p 0 ) = tee+ £GT v Ge (2.13) 

= J d 3 k £(k)G(p 0 ,k)£ T (k) + J d 3 k 1 d 3 k 2 f(fci)G(p 0 , ki)T n (p 0 , k 1: k 2 ) G(p 0 , k 2 )( r (k 2 ) ■ 

Thus, the final expression of the amplitude T has the form 

T = T n + E T XE. (2.14) 

with 


S(p 0 ,p) = f(l + GT n ) = £(p) + J d 3 k£(k)G(p 0 ,k)T 7T (p 0 ,k,p) . (2.15) 
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FIG. 1: Building blocks of the scattering amplitude. The first, second and third lines represent T w , 
S and £, respectively. The solid and dashed lines correspond to nucleons and pions, respectively. 
The filled circles represent £ and . 


III. RENORMALIZATION OF THE SCATTERING AMPLITUDE 

The expression for the scattering amplitude in Eq. (12.141) contains UV divergences. We 
perform renormalization by applying the BPHZ procedure, i.e. we subtract all divergences 
and sub-divergences of the loop diagrams and replace the LECs with their renormalized, 
finite values. In general, in renormalizable theories, subtractive renormalization can be 
realized by counter terms in the Lagrangian. Chiral effective field theory is renormalizable 
in the sense of effective field theories, i.e. all divergences can be absorbed into redefinition of 
an infinite number of counter terms. To realize subtractive renormalization in the considered 
problem, we would need to include the contributions of an infinite number of counter terms 
of the effective Lagrangian. Although this is possible for the case at hand by considering 
energy-dependent counter terms, here we only show explicitly one momentum- and energy- 
independent counter term Sz and write the contact interaction potential in a separable form 

V c (p',p) =C + C 2 (p' 2 + p 2 ) = ( 1, p’ 2 + Sz ) ( ^ ° 2 ) ( p2 | §z ) • (3.1) 

The new parameter is expressed as 

C = C - 2 C 2 Sz. (3.2) 

Thus, the contact-interaction potential has the form 

£(p) = (£i(p),6(p)) T = (1, p 2 + Sz) T . (3.3) 

The various terms contributing to the amplitude T are visualized diagrammatically in 
Fig. ID in terms of the corresponding building blocks, where in the first line the amplitude T n 
is shown. The second line represents 5, while the analogous diagrams for S T are not shown 
explicitly. The third line depicts the quantity £ which contributes to A, see Eq. (12.12[) . 
To obtain the amplitude using Eq. (12.14)1 . we first perform subtractive renormalization and 
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afterwards calculate numerically the remaining finite expressions for the quantities T, r , S, 
and X. In the following, we describe in detail how these quantities are renormalized. Since 
the amplitude T n is finite (the ultraviolet regularity of the equation for T n (12.8ft follows from 
the asymptotics (I2.5jl ). we begin our discussion with the subtractive renormalization of S. 
By writing Z(p 0 ,p) as a perturbative series as shown in Fig. [H 

Z = £ + £GVt' + ZGV w GV w + -- - , (3.4) 

it is easily seen that Z(p 0 ,p) = (Zi(p 0 ,p), Z 2 (p 0 ,p)) T can be obtained by solving the integral 
equation 

S = e + SG14. (3.5) 

This expression defines a system of equations for the quantities Zi t2 (po,p), which, using the 
explicit form of £(p) from Eq. (13 .3 p . can be written as 


5i(po,p) = 1 + J ^kZ^k) G(p 0 , k) V n (k,p) , (3.6) 

Z 2 (p 0 ,p) = p 2 + Sz + J d 3 kZ 2 (k) G(p 0 ,k)V n (k,p) . (3.7) 


The equation for Si (po,p) is free of ultraviolet divergences, see Eq. (I2.5[) . and has an ul- 

traviolet behavior Zi[po,p) ss const. On the other hand, to identify the divergences in 
Z 2 (po,p), it is convenient to consider iterations of Eq. (j3.7[) 


2 (p) = p +8z+ / d k (k + 5z) G(p 0 ,k)V n (k,p) 


d 3 k d 3 l {k 2 + Sz) G(po, k) V n (k, l ) G(po, /)! G(l, p) + 


(3.8) 


Remembering the definition of G(po, k) in Eq. (12.2[) . we simplify 


k 2 G(p 0 ,k ) = q 2 G(p 0 ,k ) 


7n z p 0 + 


\fk 2 


+ m 2 


2(2 7 r) 3 2 -)_ 777,2 

= q 2 G(p 0 ,k) + G(p 0 , k). 


(3.9) 


Substituting the above expression into Eq. (13.81) and re-organizing the perturbative series 
we obtain 


-2(po,p) = P 2 + q 2 J d 3 kG(p 0 ,k)V n (k,p ) 

+ q 2 J d 3 kd 3 lG(p 0: k)V„(kJ)G(p 0: l)V n (l,p) + ••• 

+ + J d 3 kG(p 0 , k) V n (k,p) 

+ [ d 3 kd 3 l Sz + G(po,k)V n (k,l) G(p 0 ,l)V w (l,p) + 


From this equation it is easily seen that Z 2 (po,p) can be written in the form 

S 2 (Po,p) =p 2 + q 2 l^i(Po,p) - 1] + 'MPo,p), 


(3.10) 


(3.11) 
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where the quantity 'h 7r satisfies the equation 


^Apo,p) = &r(Po,iO + 


d 3 k^ 7V (p 0 ,k)G(p 0 ,k) V n (k,p ), 


(3.12) 


with 

£n(Po,p) = I d 3 kG(p 0 , k)V v {k,p) + Sz, (3.13) 

or symbolically 

= + & = £iGV ir + 8z. (3.14) 

The G 14r term in £„. contains logarithmic divergence, which can be removed by adjusting 
the one-loop counter-term Sz to 


6z 


d 3 k G(m , k) V n (k, 0), 


(3.15) 


so that the quantities £,,-, and S 2 become finite. Moreover from Eq. (12.51) it follows that 


^ P =°°0(\np), ^ ^O(lnp). 

We now proceed with the renormalization of T, which in our scheme reduces to a subtrac¬ 
tive renormalization of E. The term in £ G £ r (the first diagram in the right-hand side of the 
third line of Fig. [TO, which is 8z- independent, contains divergences with energy-dependent 
coefficients. These divergences can be consistently subtracted using the BPHZ prescrip¬ 
tion. Those terms in £ G £ T which contain 8z linearly cancel the sub-divergences in the 
8 ^-independent part of the two-loop diagram £ G V n G £ r contained in the £ G T^G £ T part of 
the quantity E (second diagram in the righthand side of the third line of Fig. [I]). The overall 
divergence of the two-loop diagram £ G V n G £ 7 requires an additional BPHZ subtraction. 
Terms in £G£ T which contain Sz quadratically cancel the two-loop sub-divergence in the 
8 ^-independent part of the three-loop diagram £ G V n G V n G £ T (the last explicitly shown 
diagram in the third line of Fig. [[]). In addition, all ^-dependent parts of £G£ T require 
an additional subtraction of overall divergences. All other divergences appearing in the 
loop expansion of £ GT 7r G£ T are canceled automatically by contributions of the 8z counter 
term. For example, the one-loop sub-divergences of the (^-independent part of the three 
loop diagram £ G V^G V^G £ T are canceled by those expressions generated by the diagram 
£ G 14G£ T , which are linear in 8z. In the following, we provide the explicit expressions 
needed to compute the quantity E in Eq. (12.131) and define the corresponding subtractions. 
It is convenient to split E into three terms 


p—yoo 


E = E 0 + E® nite + E 


''div 


(3.16) 


The term E 0 contains only “pionless” contributions: 


So = iGf 


_ ( h(q) h(q) 
S ”° ~ l h(q) h(q) 


(3.17) 


where we have introduced the integrals 


{Io(q),l 2 (q),h(q)}= / { 1 , k 2 , (k 2 ) 2 } G(p 0 , k) 


(3.18) 
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We subtract the infinite local (polynomial in po — m) terms from these integrals to make 
them finite, so that /,(g) (i = 0,2,4) are replaced with the subtracted //*(//, g) defined as 


I* {»,<!) 


q) 

ifi^q) 


io(q) - io(w) 

9 

m 

87r 2 ^ Jm 2 — /i 2 



1 Q 

- m 

m 


+ 7T m 


— 7T m 


h{q) ~ q 2 hM - J d 3 k G(p 0} k ) = q 2 J^(/i, g) , 

^(g) - - J d 3 k ( k 2 + g 2 ) G(po, A:) = g 4 /*(//, g). 


(3.19) 


The subtracted integrals depend on the choice of the subtraction point /i. In principle, one 
has an additional freedom in fixing finite terms polynomial in p 0 — m. However, it leads to 
higher-order effects, therefore we will only study the p dependence of obtained results. 

The remaining terms of E are split into the finite and divergent parts, E® mte and Sj v , 
which are given by 


E 


finite 

7T 


(A £)(iGV r GS,+ f° jM&GUG* 

(0 K, 


(3.20) 


and 

Ef = (J 2 M£iG&+ A A {&GKG6 + 2 fe6G6+e,G^}. (3.21) 


It is straightforward to show using Eqs. (13. 91) . (13.111) and (j3 .1 4j) that two expressions for E 
given by Eq. (12.131) and Eq. (I3.16p are identical. Note that the divergent part Ej v contains 
only a finite number of iterations of the OPE potential (up to three loops as shown in Fig. [[]). 
All nonperturbative effects due to OPE are included in E® mte . 

Again, following the BPHZ procedure, we subtract the infinite local terms containing 
overall divergencies from E)( lv of the following form 


<5E, = (J 2 1 ?2 )fiG ! o6r+ (0 1 ) + 2<fc£ 1 G£ 1 + 6 r G 0 & r }, (3.22) 

where Go(k) = G(p 0 = m, k), so that the full subtracted result for E reads 


E R = (A ’4) {/„ fi (ft9)+6Gl4GS 1 }+ 0 2 M {(i(G-G„K, + £ 1 Gy,Gf,} 

+ ( ) '| ) {«.- (O - C 0 ) + C.. C V, G '!>,} . (3.23) 

The finiteness of E fi can be shown using ultraviolet behavior of Hi, ^ and 40 considered 
above. Performing subtractions in the spirit of chiral effective field theory we were supposed 
to also expand in powers of the pion mass, which we have not done here. However, the 
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non-analytic dependence of the resulting subtraction terms on the pion mass is of a higher 
order relative to the accuracy of our calculation. Note also that our final perturbative result 
which will be discussed in section IIVI depends on the choice of the renormalization scheme. 
This dependence is also of higher order. 

In final finite expressions for X we substitute the finite renormalized couplings for C and 

C 2 3 : 


X R 


(C R ) 1 - E n 


(3.24) 


Note that X R is not equal to X, because not all of the divergencies can be absorbed by 
means of redefinition of two available low-energy constants. We parameterize effectively the 
dependence of the result on the renormalization scheme by exploiting the freedom to choose 
the subtraction point /i. 


IV. THE CHOICE OF RENORMALIZATION CONDITIONS AND NUMERICAL 
RESULTS 

We are now in the position to specify the choice of renormalization conditions which, 
in the case at hand, translates into specifying the subtraction point /i. Notice that we 
have already made a specific choice for the subtractions of the integrals J 2 (g) and / 4 (g) in 
Eq. (I3.18p . It is useful to recall the key aspects of renormalization in the simple case of 
pionless EFT corresponding to = 0, see e.g. Refs. .10, 42], before dealing with the more 
complicated pionfull approach. To be specific, consider the NN S-wave scattering amplitude 
corresponding to the contact interaction potential of Eq. (12. lOj) . 


T CO nt — 


2 {C'fm 2 [1(g) g 4 - 2 g 2 / 2 (g) + J 4 (g)] + 4C 2 g 2 + 2 C} 
Clm 4 [/(g)/ 4 (g) - J 2 (g) 2 ] + 4C 2 m 2 / 2 (g) + 2 CI{q)w? - 4 


(4.1) 


To make the following discussion more transparent, we restrict ourselves to the leading 
nonrelativistic approximation so that the above expression takes the form 


T 
cont 


where 


Cfm [g 4 J(g) - 2 g 2 J 2 (g) + J 4 (g)] + 2 C 2 g 2 + C 
mJ(q ) [C%mJ 4 (q) + C] - [C 2 mJ 2 {q ) - 1] 2 

T Tf \ T / \ T f u 2 [ d 3 k {1 ,fc 2 , (fc 2 ) 2 } 


(4.2) 


(4.3) 


We have verified via explicit calculations that the omitted 1/m-correct ions are heavily sup¬ 
pressed for the problem at hand and of no relevance for the forthcoming discussion. Using 
dimensional regularization we express J 2 (g 2 ) and J 4 (g 2 ) in terms of J(g 2 ) as 


Mq 2 ) = q 2 J(q 2 ), Mq 2 ) = q 4 J(q 2 ) 

and subtract J(g 2 ) at g 2 = —p 2 obtaining 


j R (r) = 


iq + p 
2nm 


(4.4) 


(4.5) 


It is not possible to disentangle the D term from the fitted value of C R . 
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FIG. 2: Renormalization scale dependence of the phase shift in 1 So partial wave NN scattering 
emerging by non-perturbative inclusion of the NLO contact interaction in pionless EFT. Left panel 
corresponds to a natural scattering length while right panel shows the case of an unnaturally large 
scattering length. Circles (color online: red) on both panels refer to the synthetic data as described 
in the text while the dashed curves with increasing dash length correspond to /i = 1, 50,100, 200,400 
MeV. 


Subtracted integrals J B (q 2 ) and J B (q 2 ) are obtained by replacing J(q 2 ) in Eq. il 1 . Hi by 
J B (q 2 ) specified in Eq. (14.51) . The renormalized amplitude is then given by Eq. (j4.2[) with 
the divergent integrals being replaced by their subtracted values and the bare LECs C and 
C ‘2 being replaced by the renormalized ones C B (n) and C B (/i) (HI- Using the scattering 
length a and effective range r to determine these two LECs, we obtain the renormalized 
expression for the effective range function q cot 8 in terms of observable quantities 


q cot 8 


— a 2 rfi q 2 + 2 a/i — 2 
a ( arq 2 — 2 a/r + 2 ) 


i 2 2 4 

1 rq ar^q 

—|—— H---b 

a 2 A{a/i — 1) 


a 2 r 3 q e 


8(a/i — l) 2 


+ .... 


(4.6) 


Notice that the resulting expression is explicitly yU-dependent. This is because the UV 
divergencies emerging from iterations of the LS equation require counter terms beyond the 
truncated potential unless the C 2 and higher-order interactions are treated in perturbation 
theory. For the natural case describing a perturbative scenario corresponding to a ~ A -1 , 
r ~ A -1 , ..., with A being the hard scale of the order of A ~ M n , it is appropriate to 
choose the subtraction scale /i of the order of the soft scale in the problem, i.e. of the order 
of external momenta of the nucleons /i ~ q -C A. This ensures that the values of the shape 
parameters v % in Eq. (14.61) scale with the corresponding powers of A, so that the residual 
/i-dependence in the amplitude is beyond the accuracy of the NLO approximation. This is 
visualized in the left panel of Fig. [21 where the NLO pionless EFT predictions for the phase 
shift for the case of a = — M~ x fm and r = M~ l fin are shown as a function of laboratory 
energy E lab for different choices of the subtraction scale /i. The ’’data” in Fig. [2] correspond 
to the effective range approximation with all shape coefficients set to zero. Notice that 
choosing /j, of the order of the hard scale also results in a valid low-energy expansion of the 
effective-range function as visualized in the figure. 

On the other hand, for the unnatural case describing the non-perturbative situation of 
a system being close to the unitary limit and corresponding to very large values of the 


10 















scattering length, a —> oo, one obtains from Eq. (14.61) the expansion 


q cot 5 


rji q 2 
rq 2 — 2/i 


2 2 4 S 

r_^_ r_^_ 

2 4/i 8/i 2 


(4.7) 


As it is clear from Eq. (14.71) . the generated (scheme-dependent) coefficients of the effective 
range expansion will be unnaturally large if one chooses /i of the order of the soft scale in the 
problem. In the KSW approach, one compensates for these large contributions by taking 
them into account perturbatively and canceling against the contributions of the correspond¬ 
ing higher-order contact interactions which are also assumed to be unnaturally large (i.e. the 
assumed scaling of the corresponding LECs involves powers of the soft scale). We solve the 
problem of the unnaturally large scattering length by choosing the subtraction point /a of the 
order of the hard scale in the problem, /i ~ A. This guarantees that no large contributions in 
the induced coefficients of the effective range expansion are generated and the ^-dependence 
of the scattering amplitude is indeed beyond the order one is working at, see Eq. (14.61) . This 
is visualized in the right panel of Fig. [21 Notice that choosing /j, -C A leads to strong distor¬ 
tions in the phase shifts and thus considerably restricts the range of applicability of pionless 
EFT which is expected to be valid for energies up to Ej a b = M 2 /(2m) ~ 10.5 MeV. Before 
turning to the pionfull EFT we are actually interested in, it is important to emphasize that 
the subtraction scale /i should also not be chosen to be significantly larger than the corre¬ 
sponding hard scale in the problem in order to keep /x-dependent terms beyond the accuracy 
of the calculation. This feature cannot be illustrated in the considered example of pionless 
EFT, where taking the limit /a —» oo simply leads to vanishing shape parameters. In the 
presence of a long-range interaction, the induced /i-dependent contributions in Eq. (14.6j) will, 
in general, involve the mass scale associated with the long-range interaction and positive 
powers of /a. Choosing /jl^> A will then enhance the scheme-dependent contributions, which 
are nominally of a higher order, and spoil the predictive power of a theory. An explicit 
example of such a “peratization” is considered in Ref. [10:]. 

After these introductory remarks, we are now in the position to present our results for the 
phase shift at NLO in chiral EFT. We employ the exact isospin symmetry as appropriate 
at LO and use the following values for the LECs entering the OPE potential 


M n = 138 MeV, F n = 92.4 MeV, g A = 1.267. (4.8) 

The numerical value of the renormalized LECs C r {/j) and C R (n) are determined from a fit 
to the neutron-proton 1 Sq phase shift of the Nijmegen partial-wave analysis (PWA) [43| in 
the energy range of 0... 50 MeV for several choices of the subtraction point /x as discussed 
above. The resulting phase shifts for different choices of /x are plotted in Fig. [3j As expected 
and explained at the beginning of this section, we do observe some residual /i-dependence 
of the predicted phase shifts which gets strongly enhanced if one chooses /x of the order or 
smaller than M n . On the other hand, for the appropriate choice of the subtraction scale 
/i ~ A, where the hard scale A can be realistically estimated to be of the order of A ~ 600-700 
MeV 4 , the dependence on /a appears to be moderate, and the predicted energy dependence 


4 These values are of the order of the masses of the sigma and rho mesons which phenomenolo gica lly are 
known to yield the most important short-range contribution to the nucleon-nucleon potential 44L This 

A 


estimation also agrees well with the findings of chiral EFT calculations utilizing a finite cutoff B 
deeper discussion on the breakdown scale of nuclear chiral EFT can be found in Ref. 17]. 
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FIG. 3: Neutron-proton 1 5o phase shifts versus the energy in the laboratory frame. Circles (red) 
correspond to the Nijmegen PWA [43]. Dotted line represents the LO result. Curves with different 
dashing (colors) correspond to non-perturbative inclusion of the NLO contact interaction potential 
for fi = 50,100, 300, 500, 700, 850 and 900 MeV respectively. 


of the phase shift shows a good agreement with the Nijmegen PWA. Choosing /i ~ A, the 
observed /i-dependence of the phase shift is considerably smaller than the difference between 
the LO and NLO results and can serve as an estimation of the size of corrections beyond 
NLO, i.e. it defines the lower bound for the theoretical uncertainty of our calculation. It is 
especially comforting to see that the spread of predictions for the subtraction point chosen 
in the range of /i = 500... 900 MeV matches very well the estimated theoretical accuracy 
at NLO in calculations based on a finite cutoff, see Fig. 9 of Ref. [38j 

It is also interesting to address the question of perturbativeness of the subleading short- 
range interaction within our scheme. Given that calculating phase shifts always relies on 
some kind of unitarization procedure, it is more appropriate to address this issue by looking 
at the scattering amplitude directly. To be specific, consider the ratio R(E\ a b) defined as 


R{E lab) 


r NLO (£ lab ) 

|T LO (A lab )| 


(4.9) 


where T LO and T NLO denote the T-matrix calculated at LO and up to NLO, respectively. In 
Fig. [4] we show by the solid line the quantity R based on the nonperturbative inclusion of the 
subleading contact interaction as described above and corresponding to the choice of /i = 850 
MeV . 5 Notice that the subleading contribution to the amplitude becomes comparable in 
size with the leading one at higher energies, the feature that could have been expected by 
looking at the LO prediction for the 1 So phase shift. We also plot in this figure the ratio R 
resulting from the inclusion of the subleading contact interaction in first-order perturbation 


5 Note that the ratio R is not equal to 1 at threshold, because in our subtraction scheme diagrams containing 
NLO contact interactions but no overall divergencies are not subtracted. 
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FIG. 4: The ratio R defined in Eq. (14.91) versus the energy in the laboratory frame. The solid 
(dashed) line shows the result based on T NLO corresponding to the nonperturbative (perturbative) 
inclusion of the subleading contact interaction. In both cases, the subtraction scale is set to /r = 850 
MeV. 


TABLE I: Predictions for the coefficients in the effective range expansion of the 1 So phase shifts 
(low-energy theorems) at LO and NLO in the modified Weinberg approach in comparison with 


the NLO KSW predictions of Ref. 14] and empirical numbers extracted from the Nijmegen PWA 
45l . 46], For the NLO Weinberg results, we show the predictions corresponding to the variation of 
the subtraction point in the range of fi = 500 ... 900 MeV. The errors quoted for the LO predictions 
refer to the uncertainty in the numerical extraction of the coefficients [l|. 



a [fm] 

r [fm] 

V 2 [fill 3 ] 

V 3 [fm 5 ] 

V 4 [fm 7 ] 

LO, Ref. [1] 

fit 

1.50 

-1.9 

8 . 6 ( 8 ) 

-37(10) 

NLO, nonperturb. C 2 

fit 

fit 

-0.55... -0.61 

5.1... 5.5 

-29.6... - 30.8 

NLO, perturbative C 2 

fit 

fit 

-0.51... -0.57 

4.5... 4.7 

-28.8... -29.8 

NLO KSW, Ref. [14] 

fit 

fit 

-3.3 

18 

-108 

Nijmegen PWA 

-23.7 

2.67 

-0.5 

4.0 

-20 


theory for the same choice of /i (and using the same values for the renormalized low-energy 
constants as determined in the nonperturbative calculation). This shows clearly that it is 
advantageous to include the subleading contact interaction nonperturbatively within the 
employed framework for energies of about E\^ ~ 50 MeV and higher. 

Last but not least, we also give the predictions for the coefficients in the effective range 


expansion which may be regarded as low-energy theorems (LETs), see 10, LL4] for more 
details. In table HI the LETs in the KSW and Weinberg approaches are confronted with 
the results of the Nijmegen PWA for the 1 S'o partial wave, respectively. We observe a clear 
improvement in the reproduction of the LETs when going from LO to NLO. Notice that 
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the nonperturbative treatment of the subleading contact interaction appears to have minor 
effect for the LETs. It should, however, be emphasized that the extraction of the coefficients 
in the effective range expansion requires performing a unitarization of the amplitude which 
provides a partial resummation of (^-contributions. 


V. SUMMARY AND CONCLUSIONS 


In this paper we have considered nucleon-nucleon scattering in the 1 S'o partial wave 
within the modified Weinberg approach. The integral equation based on the leading-order 
potential, which consists of the momentum- and energy-independent contact interaction 
and the OPE potential, is renormalizable and was studied in Ref. [lj. The observed large 
discrepancy between the LO EFT results and the 1 So phase shift of the Nijmegen PWA, 
which starts already at rather low energies, indicates that at least some parts of the higher- 
order contributions to the effective potential need to be included nonperturbatively. Here we 
assumed that only the short range part of the NLO potential (in standard Weinberg power 
counting) needs to be treated non-perturbatively. It involves only the contact interaction 
terms quadratic in momenta and the pion mass. This makes it possible to perform the 
subtractive renormalization explicitly in non-perturbative expressions. 

The pertinent results of our study can be summarized as follows: 


• We have carried out subtractive renormalization of the scattering amplitude based on 
the potential involving OPE as well as the leading and subleading contact interactions 
using the framework of Ref. jlj and without relying on perturbation theory. 

• The resulting renormalized integral equations for the scattering amplitude have been 
solved numerically and the values of the renormalized low-energy constants C R (n) and 
C' R (n) were determined from a fit to phase shifts of the Nijmegen PWA for different 
choices of the subtraction point /i. 


• We discussed the issue of the proper choice of renormalization conditions in our scheme 
and have argued that the observed large value of the scattering length requires choosing 
the scale / 1 , which corresponds to the renormalization of the LO contact interaction, 
of the order of the hard scale in the problem. 


The resulting preditions for the energy dependence of the 1 So phase shift are in a good 
agreement with the Nijmegen PWA. Moreover, the observed dependence of the phase 
shifts on the subtraction point /a chosen in the range of /1 = 500... 900 MeV ag rees 
well with the theoretical accuracy at NLO estimated in calculations of Refs. [1 381 39l 
based on the standard non-relativistic framework with a finite cutoff. 


• We have also addressed perturbativeness of the subleading contact interaction within 
our scheme. We found that it is advantageous to treat the subleading contact inter¬ 
action nonperturbatively at energies of about E^b ~ 50 MeV and higher. 

• Finally, we have looked at the low-energy theorems for the coefficients in the effective 
range expansion and found a clear improvement when going from LO to NLO. 


The results of our work open the way to perform higher-order calculations within the mod¬ 
ified Weinberg approach proposed in Ref. [lj. As a next step, the role of the two-pion 
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exchange potential needs to be investigated and the extension to other partial waves has to 
be performed. Work along these lines is in progress. 
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